Auger recombination and carrier multiplication in embedded silicon and germanium 

nanocrystals 
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For Si and Ge nanocrystals (NCs) embedded in wide band-gap matrices, Auger recombination 
(AR) and carrier multiplication (CM) lifetimes are computed exactly in a three-dimensional real 
space grid using empirical pseudopotential wave functions. Our results in support of recent experi- 
mental data offer new predictions. We extract simple Auger constants valid for NCs. We show that 
both Si and Ge NCs can benefit from photovoltaic efficiency improvement via CM due to the fact 
that under an optical excitation exceeding twice the band gap energy, the electrons gain lion's share 
from the total excess energy and can cause a CM. We predict that CM becomes especially efficient 
for hot electrons with an excess energy of about 1 eV above the CM threshold. 

PACS numbers: 72.20.Jv, 73.22.-f 



I. INTRODUCTION 

The Group IV semiconductor, Si and to a lesser ex- 
tend, Ge have been indispensable for the electronics and 
photovoltaics industry. The recent research efforts have 
shown that their nanocrystals (NCs) bring new features 
which fortify their stands. For instance, NCs can turn 
these indirect band-gap bulk materials into light emit- 
ter si or offer increased efficiencies in solar cells 3 The lat- 
ter has been demonstrated in a very recent experimental 
study^ by significantly increasing the solar cell efficiency 
in colloidal Si NCs due to carrier multiplication (CM) 
which enables multiple exciton generation in response to 
a single absorbed photon4£ Similarly, the inverse pro- 
cess, Auger recombination (AR) is also operational and 
it introduces a competing mechanism to CM which can 
potentially diminish the solar cell efficiency and in the 
case of light sources it degrades the performance by in- 
flating the nonradiative carrier relaxation rate£ 

The utilization and full control of both CM and AR 
require a rigorous theoretical understanding. The pi- 
oneering series of theory publications on the AR in Si 
NCs belong to a single group based on an atomistic tight 
binding approach*^^ Unfortunately, they only consid- 
ered hydrogen passivated Si NCs without addressing the 
shape effects. Moreover, their results do not reveal a size- 
scaling for AR but rather a scattered behavior over a wide 
band of lifetimes in the range from few picoseconds to few 
nanoseconds as the NC diameter changes from 2 to 4 nm. 
In the past decade, no further theoretical assessment of 
AR in Si NCs was put forward. In this context, the Ge 
NCs have not received any attention although with their 
narrower effective band gap, they can benefit more from 
the low-energy part of the solar spectrum in conjunction 
with CM for increasing the efficiency. 

In this work we provide a theoretical account of AR 
and CM in Si and Ge NCs which reveals their size, shape 
and energy dependence. Another important feature of 
this work, unlike commonly studied hydrogen-passivated 
NCs is that we consider NCs embedded in a wide band- 
gap matrix which is essential for the solid-state device 




FIG. 1: (Color online) AR in NCs: (a) excited electron, 
(b) excited hole; the solid and dashed arrows refer to direct 
and exchange processes. Biexciton type AR and its inverse 
process CM: (c) excited electron (d) excited hole. 



realizations. Similar to the classification of Wang et al. in 
their theoretical work on Coulombic excitations in CdSe 
NCsjiS we consider different possibilities of AR as shown 
in Fig. [TJ We use the type of the exicted carrier as the 
discriminating label, hence we have the excited electron 
(Fig- Hl^a)) and the excited hole (Fig.QJb)) AR and their 
biexciton variants (Fig. [TJc) and (d)). All of these have 
corresponding CM processes taking place in the reverse 
direction but only the CM in Fig. QJc) and (d) are studied 
as they can be optically induced. 



II. THEORETICAL DETAILS 

Both AR and CM require an accurate electronic struc- 
ture over a wide energy band extending up to at least 
3-4 eV below (above) the highest occupied molecu- 
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FIG. 2: (Color online) (a) Embedded NC in a supercell, core 
atoms of a (b) spherical, (c) oblate and (d) prolate ellipsoidal 
NC. 



lar orbital-HOMO (lowest unoccupied molecular orbital- 
LUMO). Another constraint is to incorporate several 
thousands of core and host matrix atoms within a su- 
percell (see Fig. [5] (a)). To meet these requirements we 
have employed the linear combination of bulk bands basis 
within the empirical pseudopotential framework which 
can handle thousands-of-atom systems both with suffi- 
cient accuracy and efficiency over a large energy win- 
dow Details regarding its performance and the imple- 
mentation such as the wide band-gap host matrix can be 
found in Ref. [ljj We should mention that Califano et 
al. have successfully employed a very similar theoretical 
approach in order to explain the hole relaxation in CdSe 
NCs.i 3 - 

After solving the atomistic empirical pseudopotential 
Hamiltonian for the energy levels and the wave functions, 
the AR and CM probability can be extracted using the 
Fermi's golden rule, 
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where ipi and ipf are respective initial and final config- 
urations with the corresponding energies Ei and Ef, re- 
spectively, and r is the level broadening parameter which 
is taken as 10 meV. However sensitivity to this parame- 
ter is also considered in this work. The spin-conserving 
screened Coulomb potential is given by V c (i"i,r 2 ) = 
e 2 /e(rx, r 2 )|ri — r 2 |. The subject of screened Coulomb 
interaction in NCs is an active source of debate; re- 
cent publications predict reduced screenin g 14 i 15 whereas, 
other theoretical investigations^^ have concluded that 
the inverse dielectric function is bulklike inside the NC. 
Therefore, we followiS and use 
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ro(ri)m(r 2 ), (2) 



as the inverse dielectric function, where, the so-called 
mask function m(r) is set to 1 when r inside of the NC 
and when r outside of the NC. 

Expressing the initial and final states of the AR shown 
in Fig. Q] (a) or (b) by using the Slater determinant, the 
matrix elements ((ipi \V c (\ri, r 2 |)| ipf)) can be calculated 



M(i,j;k,l) = ±5 J J <P* k (ri)<Pt(r2)V c (r 1 ,v 2 ) 

x(^(n)^(r 2 ) - <Mr 2 )<Mn))d 3 rid 3 e3) 

here the labels i, j and k and I refer respectively to the 
initial and final states which also include the spin and V 
is volume of the supercell. 
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FIG. 3: (Color online) AR lifetimes for (a) excited electron, 
(b) excited hole, and (c) biexciton types in Si NCs, and (d) ex- 
cited electron, (e) excited hole, and (f) biexciton types in Ge 
NCs. Square symbols represent AR lifetimes with interface 
smearing, and dashed lines show AR lifetimes calculated from 
our proposed C values. Spherical symbols in (a) and (b) rep- 
resent AR lifetimes in Si NCs with abrupt interfaces. 

These matrix elements, M(i,j;k,l) are computed ex- 
actly in a three-dimensional real space grid without re- 
sorting to any envelope approximation. The number of 
final states are determined setting the final state win- 
dow to ±7r around the exact conserved energy E^(= 
Ej + Ei — Ei). For the initial states i and I, Boltzmann 
average is taken into account due to thermal excitations. 
The same should apply to the other initial state j, how- 
ever, as a safe but computationally very rewarding sim- 
plification it is kept fixed at LUMO for the excited elec- 
tron (EE), and at HOMO for the excited hole (EH) type 
AR. 



III. RESULTS AND DISCUSSIONS 

We first apply this formalism to spherical NCs (see 
Fig. [2b)) having abrupt interfaces. The corresponding 
AR lifetimes for EE and EH processes are plotted as a 
function of NC diameter in Fig. [3] (a) and (b). The C% v 
point symmetry of the NCs in the case of abrupt interface 
between NC core and the matrix causes oscillations in the 
physical quantities such as the state splittings and the 
density of states with respect to NC diameter^ When 
we account for the interface transition region between the 



NC and host matrix^ we observe that these strong oscil- 
lations in the size dependence of AR are highly reduced 
for Si and Ge NCs (cf. Fig. [3]). The interface region 
especially affects the excited state wave functions and 
the final state density of states and it makes our model 
more realistic for both Si and Ge NCs. As an observa- 
tion of practical importance, we can reproduce our data 
remarkably well (cf. Fig. [3J using the simple expression 
1/t = Cn 2 , with an Auger coefficient of C = 1 x 10~ 30 
cnrV 1 for Si NCs and C = 1.5 X 1CT 30 cm 6 s~ 1 for Ge 
NCs, where n = 2/V NC is the carrier density within the 
NC such that there should be two electrons or holes to 
initiate an AR. 

The other two types of AR shown in Fig. [3] (c) and (f ) 
refer to biexciton recombinations. This process becomes 
particularly important under high carrier densities such 
as in NC lasers or in solar cells under concentrated sun- 
light. Its probability can be expressed in terms of EE 
and EH type AR as^ 1/ Txx = 2/r EE + 2/r EH where r EB 
and r EH are EE and EH lifetimes. Fig. [3] (c) and (f) 
compares the computed biexciton type AR for Si and 
Ge NCs with the expression 1/t = Cn 2 where the value 
C = 4xl0~ 30 cm 6 s~ 1 and 6x 10 _30 cm 6 s _1 are used which 
are obtained from the previous C values extracted for EE 
and EH processes together with the r xx expression. For 
Si NC case, our calculated value at 3 nm diameter agrees 
reasonably well with the experimental photoluminescense 
decay time of about 105 ps which was attributed to AR. 18 
In Fig. [3J (c) and (f), we also demonstrate the fact that 
a choice of V = 5 meV does not introduce any marked 
deviation from the case of T = 10 meV as used in this 
work for both Si and Ge NCs. This parameter test auto- 
matically checks the sensitivity to the final state energy 
window chosen as ±7r. 



TABLE I: AR lifetimes for different ellipsoidal shapes of Si 
NCs with diameters of 1.63 and 2.16 nm. 





Spherical 


Prolate 


Oblate 


D (nm) 


1.63 2.16 


1.63 2.16 


1.63 2.16 


EE (ps) 


40 541 


32 103 


36 121 


EH (ps) 


267 430 


74 76 


26 139 



Next, we demonstrate the effects of deviation from 
sphericity on Si NCs. We consider both oblate (Fig.^c)) 
and prolate (Fig. [DJd)) ellipsoidal Si NCs described by 
the ellipticity value of e=0.85. For the comparison pur- 
poses, we preserve the same number of atoms used in 
spherical NCs of diameters 1.63 and 2.16 nm. The re- 
sults listed in Table U indicate that the spherical NC has 
a lower Auger rate than the aspherical shapes. This can 
be reconciled as follows: in the case of either prolate or 
oblate NC, the electronic structure is modified in such 
a way that the number of final states is increased, fur- 
thermore, a coalescence of the states around the HOMO 
and LUMO occurs. A similar effect was also observed in 
the asphericity-induced enhancement of Auger thermal- 
ization of holes in CdSe NCs.- 3 However, we should note 



that the shape effects are not pronounced. 

In their work on the CM in PbSe NCs, Allan and 
Delerue have deduced that such Coulombic interactions 
are primarily governed by the state-density function. 19 
Even though we agree on the importance of the density 
of states, we would like to emphasize the significant role 
of the Coulomb matrix elements. We illustrate our point 
by Fig. 31 where the average matrix element for Si and 
Ge NCs are shown. The strong size dependence leading 
to a variation over several orders of magnitude indicates 
their nontrivial role. 
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FIG. 4: (Color online) Average Coulomb matrix elements for 
(a) Si and (b) Ge NCs for EE type AR (red squares) and for 
EH type AR (black spheres). 



Regarding the CM, we first consider the inverse Auger 
process (cf., Fig. [1] (c) and (d)) for different diameters 
of the Si and Ge NCs. Therefore, we consider the im- 
pacting electron (hole) to have an energy of E gap = 
-Elumo — -Ehomo above (below) the conduction (valence) 
band edge i.e., just at the threshold energy to initiate a 
CM event. As seen in Fig. [5j EE and EH type CM life- 
times for Si and Ge NCs decrease from the few ns to 
few ps as the NC diameter decreases. However, for EE 
(EH) type CM, the small number of final states at the 
bottom of the CB (top of the VB) cause a nonmonotonic 
dependence of CM on size of the NC. 
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FIG. 5: (Color online) CM Lifetime results for (a) EE and 

(b) EH types in Si NCs embedded in SiCb and AI2O3, and 

(c) EE and (d) EH types in Ge NCs embedded in AI2O3. 

From the practical point of view, the investigation of 
the effect of excess energy on the CM under an optical ex- 
citation above the effective gap, E gap is even more impor- 
tant. Here, the correct placement of excited electron and 
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FIG. 6: (Color online) (a) Electron and hole excess energy vs 
optical pump (excitation) energy for 4 nm Si and 3 nm Ge 
NCs, (b) CM lifetime vs electron excess energy for different 
diameter of Si and Ge NCs. The horizontal arrows provide 
the 1 and 0.5 eV marks. 



tion of CM in photovoltaic applications^ Our threshold 
value of 2.8 E gap for Si NCs is somewhat higher than the 
recent experimental data£ In Fig. ® (b) we show the cor- 
responding electron-initiated CM lifetimes as a function 
of excess energy. It can be observed that CM is enhanced 
by more than two orders of magnitude within an excess 
energy of 1 eV beyond the CM threshold reaching a life- 
time of few picoseconds. 



IV. CONCLUSIONS 

To summarize, we offer a theoretical assessment of the 
two most important Coulombic excitations, AR and CM, 
in Si and Ge NCs. The Auger coefficients that we ex- 
tracted can serve for the practical needs in the utilization 
of this process. For the efficiency enhancement via CM 
in Si and Ge NCs, the prospects look positive as the hot 
electrons receive most of the excess energy and they can 
undergo a CM within few picoseconds. 



hole after optical excitation is critical. We assign the ex- 
cited electron and hole to their final states by accounting 
for all interband transitions with the given energy differ- 
ence as weighted by the dipole oscillator strength of each 
transition which is a direct measure of the probability of 
that particular event. In Fig. [5] (a) we observe that the 
electrons receive the lion's share of the total excess energy 
which is the desired case for the high efficiency utiliza- 
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